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Abstract. The problem of anomalous scaling in the model of a transverse vector field 
9i(t,x) passively advected by the non-Gaussian, correlated in time turbulent velocity 
field governed by the Navier-Stokes equation, is studied by means of the field-theoretic 
renormalization group and operator product expansion. The anomalous exponents of 
the 2n-th order structure function S^nM = x ) — 0(t,x + r)] 2n ), where 8 is the 
component of the vector field parallel to the separation r, are determined by the critical 
dimensions of the family of composite fields (operators) of the form {d9d0) 2n , which 
mix heavily in renormalization. The daunting task of the calculation of the matrices 
of their critical dimensions (whose eigenvalues determine the anomalous exponents) 
simplifies drastically in the limit of high spatial dimension, d — > oo. This allowed us 
to find the leading and correction anomalous exponents for the structure functions up 
to the order Sse- They reveal intriguing regularities, which suggest for the anomalous 
exponents simple "empiric" formulae that become practically exact for n large enough. 
Along with the explicit results for modest n, they provide the full description of the 
anomalous scaling in the model. Key words: passive vector field, turbulent advection, 
anomalous scaling, renormalization group, operator product expansion. 
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1. Introduction 

In the past two decades, much attention has been attracted by turbulent advection of 
passive scalar fields; see the review paper p] and references therein. Being of practical 
importance in itself, the problem of passive advection can be viewed as a starting point 
for studying intermittency and anomalous scaling in the fluid turbulence on the whole [2]. 
Most progress was achieved for the so-called Kraichnan's rapid-change model, in which 
the advecting velocity field Vi(x) with x = {t,x} is modelled by a Gaussian statistics 
with vanishing correlation time and prescribed correlation function (vv ) oc 5(t—t')k~ d ~^ , 
where k is the wave number, d is the dimension of space and £ is an arbitrary exponent 
with the most realistic (Kolmogorov) value £ = 4/3. The structure functions of the 
advected scalar field 9{x) in the inertial range demonstrate anomalous scaling behaviour: 

S 2n (r) = ([0(f,x) - 6(ty)} 2n ) oc r n ^ (r/L) A ", (1.1) 

that is, singular dependence on the separation r = |r| (where r = x — x') and on the 
integral turbulence scale L, characterized by an infinite set of the exponents A n . Within 
the framework of the so-called zero-mode approach, these exponents were calculated in 
the leading order of the expansions in £ [3] and 1/d [I] : 

A n = -2n(n - l)£/(d + 2) + 0(£ 2 ) = -2n(n - l)£/d + 0(l/d 2 ). (1.2) 

In [5] and subsequent papers, the field theoretic renormalization group (RG) and 
the operator product expansion (OPE) were applied to Kraichnan's model; see [6] for 
a review and the references. In that approach, the anomalous scaling emerges as a 
consequence of the existence in the corresponding OPE of certain composite fields 
("operators" in the quantum- field terminology) with negative dimensions, which are 
identified with the anomalous exponents A n . This allows one to construct a systematic 
perturbation expansion for the anomalous exponents and to calculate them up to the 
orders £ 2 [5] and £ 3 [7]. Besides the calculational efficiency, an important advantage of 
the RG+OPE approach is its relative universality: it can also be applied to the case of 
finite correlation time or non-Gaussian advecting fields. For passively advected vector 
fields, any calculation of the exponents for higher-order correlations calls for the RG 
techniques already in the 0(£) approximation. 

In this paper, we study anomalous scaling of a passive vector quantity, advected 
by a non-Gaussian velocity field, governed by the stirred Navier-Stokes (NS) equation. 
For the rapid-change velocity ensemble, a similar model was introduced and thoroughly 
studied in [8] [12]. Before explaining our motivations, which follow the same lines as 
those of refs. [8j QUI QT], let us discuss the definition of the model. 

We confine ourselves with the case of transverse (divergence-free) passive 9i(x) and 
advecting Vi(x) vector fields. Then the general advection-diffusion equation has the form 

- A(e k d k )vi + d-P = K d% + rji, V t = d t + (v k d k ), (1.3) 

where V f is the Lagrangian (Galilean covariant) derivative, V(x) is the pressure, n$ is 
the diffusivity, d 2 is the Laplace operator and rji(x) is a transverse Gaussian stirring 
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force with zero mean and covariance 

( Vi (x) Vk (x')) = 8(t - t') C lk (r/L). (1.4) 

The parameter L is an integral scale related to the stirring, and Cik is a dimensionless 
function, finite at r = and rapidly decaying for r — > oo; its precise form is unimportant. 
Due to the transversality conditions diOi = diVi = 0, the pressure can be expressed as 
the solution of the Poisson equation, 

d 2 v = {A-\)d l v k d k e l . (i.5) 

Thus the pressure term makes the dynamics (11.31) consistent with the transversality. 
The amplitude factor A in front of the "stretching term" (8kdk)vi is not fixed by the 
Galilean symmetry and thus can be arbitrary. Such general 11 A model" was introduced 
and studied in refs. [13]. The most popular special case A — 1, where the pressure term 
disappears, corresponds to magnetohydrodynamic turbulence. It was studied earlier in 
numerous papers; see e.g. refs. [H] and references therein. 

In earlier studies, the velocity field in (11.31) was described by the Kraichnan's rapid- 
change model. In this paper, we employ the stochastic NS equation: 

V t Vi = is d 2 Vi - dip + fi, (1.6) 

where Vf is the same Lagrangian derivative, p and f\ are the pressure and the transverse 
random force per unit mass. We assume for / a Gaussian distribution with zero mean 
and correlation function 

(fi(x)f j (x f )) = 5 -^p- [ rfkP 4J (k)rf / (A;)exp[ik(x-x')], (1.7) 

l Z7r J Jk>m 

where -Py(k) = 5ij — kikj/k 2 is the transverse projector, df(k) is some function of k = |k| 
and model parameters. The momentum m = 1/L, the reciprocal of the integral scale L 
related to the velocity, provides IR regularization. For simplicity, we do not distinguish 
it from the integral scale related to the scalar noise in ( 11. 4p . 

The standard RG formalism is applicable to the problem (11. 6p . (I1.7P if the 
correlation function of the random force is chosen in the power form [15] 

d f (k) = D k 4 - d - y , (1.8) 

where D > is the positive amplitude factor and the exponent < y < 4 plays the 
role of the RG expansion parameter. The most realistic value of the exponent is y — 4: 
with an appropriate choice of the amplitude, the function (I1.8P for y — >■ 4 turns to the 
delta function, df(k) ex 5(k), which corresponds to the injection of energy to the system 
owing to interaction with the largest turbulent eddies; for a more detailed justification 
see e.g. [TOl fT7] . 

In this paper we consider the model (I1.3P without the stretching term, that is, 
A = 0. Being formally a special case of the general A model, it appears exceptional in 
a few respects and requires special attention. 

The feature specific only of the A = is the symmetry with respect to the shift 
9 — t- 9 + const, because only derivatives of the field 9 enter the equation (11. 3ft . The 
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quantities of interest are the structure functions (11.11) . in which 9 should be understood 
as the component of the vector field parallel to the separation, 9 = 9iri/r: in contrast 
to ordinary correlation functions, they are also invariant with respect to the shift. As 
a consequence, all the composite operators that enter the corresponding OPE, should 
also be invariant, that is, built only of the derivatives of 9. 

For the scalar problem, the operator that determines the leading term of the inertial- 
range asymptotic behaviour (ll.ip . is unique: it has the form (di9di9) n , the n-th power 
of the local dissipation rate of the scalar field fluctuations, and its critical dimension 
gives the anomalous exponent A n is (11. 2ft ; see [5]-[7] for the detailed discussion. 

For the vector case one can construct many scalar operators of the form (d9) 2n for 
a given n, and in order to find the corresponding set of exponents and to identify the 
leading contribution, one has to consider the renormalization of the whole family, which 
implies the mixing of individual operators [51 [TP]. Renormalization of families of mixing 
composite fields and calculation of the corresponding matrices of critical dimensions 
(whose eigenvalues give the desired anomalous exponents) its rather cumbersome and 
labor-compiling task, which should be solved separately for different families (in the case 
at hand, for different n). Thus, at first sight, there is no hope to derive simple explicit 
expressions for the anomalous exponents, similar to (11.21) in the scalar case. 

In this respect, the ^4 = case of the model (jl.3p resembles the nonlinear stirred 
NS equation, where the inertial-range behavior of structure functions is believed to be 
related with the Galilean-invariant (and hence built of the velocity gradients) operators, 
which mix heavily in renormalization; see [16] and references therein. In that case, the 
full solution is not yet obtained even for the relatively simple case of the family that 
includes the square of the energy dissipation rate [18J. 

Thus the A = vector model (11.31) is of special interest: it also involves the problem 
of mixing, but now the problem is not a hopeless one: the leading terms are determined 
by finite families of composite operators, namely, those of the form (d9) 2n with all 
possible contractions of vector indices, and such families with a given n are closed with 
respect to the renormalization [HI ID]. What is more, for low values of d there are linear 
relations between the operators, which reduce drastically the number of independent 
monomials [12J. As a result, the leading anomalous exponents for d = 2 were calculated 
to the order 0(£ 2 ) for all n, and for d = 3 - to the order 0(£) for n < 9 [12]. Crucial 
simplifications also take place in the limit d — > oo [ID]: like for the scalar Kraichnan's 
model, the anomalous exponents decay as 0(l/d) for large d, so that the anomalous 
scaling disappears at d = oo; cf. (II. 2p . In order to find the leading exponents (and the 
closest corrections) it is sufficient to consider some special subset of the whole family 
(d9) 2n , and the corresponding matrix of critical dimensions can be built by a simple 
algorithm [ID]. This allowed us to find all the negative dimensions in the leading order 
of the double expansion in £ and 1/d for n as large as n = 28, which gives the anomalous 
exponents and the close corrections for the structure functions up to S 56 [HJ . All those 
results, however, referred to the Gaussian velocity field. 
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For very large n the calculations become too labor- and time compiling (mostly 
due to the diagonalization of the matrices), but they are not necessary: the large- n 
results suggest some simple empiric explicit expressions for the leading, next-to-leading, 
etc, anomalous exponents, which become practically exact as n increases. Along with 
the explicit answers for smaller n, this gives the complete description of the anomalous 
scaling in the vector model for all n and large d [IT]. 

It should be emphasized that the study of the large-ci behaviour of the fluid 
turbulence is by no means of only academic interest. It is related to the old idea of 
the expansion in 1/d, which has repeatedly been introduced in various contexts [19] 

The problem is that the ordinary perturbation theory for the stirred (stochastic) 
NS equation (that is, the perturbation expansion in the nonlinearity) is in fact an 
expansion in the Reynolds number, a parameter which tends to infinity for the fully 
developed turbulence. A similar problem is well known in the theory of critical state, 
where it is solved by means of the RG techniques; see e.g. p2]. The RG allows one 
to rearrange (to sum up) the plain perturbation series and to replace them with the 
famous e expansion, where e = 4 — d is the deviation of the spatial dimension d from 
its upper critical value d = 4. The turbulence (or, better to say, the corresponding 
stochastic models), has no upper critical dimension, and the RG expansion parameter 
has completely different meaning. As already mentioned, in Kraichnan's model its role 
is played by the exponent £, while in the RG approach to the stirred NS equation its 
analog is the exponent y in the correlator of the stirring force; see section [2j The results 
of the RG analysis of this model are reliable and internally consistent for asymptotically 
small £ or y, while the possibility of their extrapolation to the physical finite values, 
and thus their relevance for the real fluid turbulence, is sometimes called in question; 
see the discussion and the references in [23], [25] . 

One can hope that in the limit d — > oo intermittency and anomalous scaling 
disappear or acquire a simple "calculable" form and the finite-dimensional turbulence 
can be studied within the expansion around this "solvable" limit [20]. Indeed, it was 
argued (on the basis of a certain SDE-motivated ansatz for dissipative terms) that the 
Kolmogorov theory becomes exact and the multiscaling indeed disappears for d = oo 
[21~] . as it also happens for the Obukhov-Kraichnan model (THJ. What is more, for the 
latter it was possible to find the 0(l/d) contribution to the anomalous exponents [JJ; 
see the last expression in C I X . 2 1) . However, the systematic expansion in 1/d has not been 
yet constructed for that model, let alone the stochastic NS equation. 

It was suggested in |23l [24] that new progress can be achieved by combining the 
large-d limit with the RG approach and the expansion in e. In particular, it was 
noticed [TUl [23] that taking the limit d — > oo leads to serious simplifications in the 
RG calculations, especially for composite operators. In a very important paper [23J, 
scaling dimensions of all the powers of the local energy dissipation rate for the NS 
problem were calculated for d = oo to first order in y, the problem that looks unfeasible 
for finite d; see the discussion in [IS] . 
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Thus the drastic simplifications that occur in our vector model with the turbulent 
mixing provided by the NS velocity field and the simple explicit results for the anomalous 
dimensions in the leading order of the double expansion in y and 1/d, give a new strong 
support to the idea of the large-<i expansion. 

The plan of the paper is as follows. In sec. [2] we discuss the field theoretic 
formulation of our stochastic problem and its renormalization. We show that the 
corresponding RG equations have an IR attractive fixed point, which implies existence of 
IR scaling behaviour for various correlation functions. In sec. |3j inertial-range anomalous 
scaling of the structure functions is studied by means of the OPE and the role of the 
operators (dd) 2n is clarified. In sec. H]we discuss the renormalization of those operators 
and the simplifications that occur in the limit of large d. Examples are given of the 
matrices of critical dimensions for a few families of those operators. In sec. [5] the leading 
and correction anomalous exponents, determined by the eigenvalues of those matrices, 
are presented and the regularities that they reveal are discussed. These regularities 
suggest for the eigenvalues some simple "empiric" formulae that become practically 
exact for n large enough. Along with the explicit results obtained for modest n, they 
provide the full description of the anomalous scaling in the present model. Sec. |6] is 
reserved for a brief conclusion. 



2. Field theoretic formulation, renormalization and RG equations 

According to the general theorem (see e.g. [H]), the full-scale stochastic problem (11. Sp - 
i ll. 81) for A = is equivalent to the field theoretic model of the doubled set of fields 
$ = {v , v ', 6, 6'} with the action functional 

S($)=S v (V,v) + 6'D e 6'/2 + 6'{-V t + Kod 2 }6, (2.1) 

where Dq is the correlation function (11.41) of the random noise / in ( 11. 3D and S v is the 
action for the problem (jl.6j) - (jl.8p : 

S v (y\ v) = v'D v v'/2 + v' {-V t + u d 2 } v, (2.2) 

where D v is the correlation function (11.71) of the random force /j. All the integrations 
over x = {t, x} and summations over the vector indices are understood. The auxiliary 
vector fields v', 8' are also transverse, div[ = d{8i = 0, which allows to omit the pressure 
terms on the right-hand sides of expressions (I2.ip . (I2.2p . as becomes evident after the 
integration by parts. For example, 



J dt J cfx v[dip = — J dt J dx pi^div'^ = 0. 



Of course, this does not mean that the pressure contributions are unimportant: the 
fields v', 9' act as transverse projectors and select the transverse parts of the expressions 
to which they are contracted. 

The role of the coupling constants is played by the two parameters g = Dq/vq and 
Uo = kq/i>o, the analog of the inverse Prandtl number in the scalar case. By dimension, 

g oc A v and u oc A , (2.3) 
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where A is the characteristic ultraviolet (UV) momentum scale. Thus the model (12. ip . 
(12. 2 p becomes logarithmic (both the coupling constants become dimensionless) at y — 0, 
and the UV divergences manifest themselves as poles in y. 

The renormalization and RG analysis of the model (12. ip . (12. 2 p are similar to that 
of the scalar advection by the NS velocity field [261 EZ]j an d here we discuss them 
only briefly. Dimensional analysis, augmented with symmetry considerations (Galilean 
invariance and the symmetry with respect to the shift 9 — > 9 + const) shows that 
superficial UV divergences are present only in the 1-irreducible Green functions (v'v) 
and (9'9), and the corresponding counterterms have the forms v'd 2 v and 9'd 2 9. They 
can be reproduced by multiplicative renormalization of the parameters 

u = vZ v , k = kZ k , g = gn v Z g , Z g = Z~ 3 ; (2.4) 

no renormalization of the fields $ and the IR scale m is needed. Here u, g, k are 
renormalized analogs of the bare parameters uo, go, k,q and the reference scale fi is an 
additional parameter of the renormalized theory. The last relation in (12 .4p follows from 
the absence of renormalization of the amplitude Do = go^o = 9H Vl/3 in the first term 
of the action (12. 2p . The renormalization constants Zj = Zi(g,u,d,y) absorb all the UV 
divergences, so that the Green functions are UV finite (that is, finite at y — 0) when 
expressed in renormalized parameters. 
The one-loop calculation gives: 

Z'^-^m^^ht^ ^- (2 - 5) 

where S d = S d /(2n) d and S d = 2u d / 2 /T(d/2) is the surface area of the unit sphere in 
d- dimensional space. Of course, due to the passivity of the field 9, the constant Z u is 
the same as in the model (12. 2p . 

Since the model is multiplicatively renormalizable, the RG equations can be derived 
in the standard manner; see e.g. [32]. The RG equation for a certain renormalized Green 
function G R = ($...$) has the form 

{P M - lv V v + I3 g d g + Mu] G R = 0. (2.6) 

Here T> s = sd s for any variable s, u = k/v and the RG functions (the f3 functions and 
the anomalous dimensions 7) are defined as 

1f = V^\yiZ f (2.7) 

for any quantity F and 

Po = 'D li g = g[-y + &y„], f3 u = V„u = u[<y„ - ~/ K ], (2.8) 

where T>^ is the operation at fixed bare parameters and the second relations in (12. 8p 
follow from the definitions and the relations (12 .4p . It remains to note that the differential 
operator in (12. 6p is nothing but T>^ expressed in renormalized variables. 
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From (12. 5p one obtains the following explicit one-loop expressions for the anomalous 
dimensions: 

It is well known that IR asymptotic behaviour of a multiplicatively renormalizable 
field theory is governed by IR attractive fixed points of the corresponding RG equations. 
Their coordinates are found from the requirement that all the /9 functions vanish; in our 
case, (3 g = j3 u = 0. From the explicit expressions (I2.8p . (I2.9p for g it follows that the 
model (12. 2 p has a nontrivial fixed point 

9 * Bd = y W^T) +0{y2) (2 - 10) 

which is positive and IR attractive (d g (3 g > 0) for y > (of course, this fact is well 
known, see e.g. [HJ EE]). Then from the equation (3 U = and the explicit expressions 
flUD, (EU) and (I2A0D one obtains 

u *( u * + 1 ) = %zr^+°(y)- (2- 11 ) 

We are interested in the positive root of the equation (12. lip which exists for d 2 > 3. It 
is easily checked that this fixed point is IR attractive (d u (3 g = 0, d u (3 u > 0). 

Existence of an IR attractive fixed point in the physical range of parameters (w* > 0, 
g* > 0) means that the Green functions of the full model (12. ip . (12.21) exhibit self-similar 
(scaling) behaviour in the IR asymptotic range. The corresponding critical dimensions 
A[F] = Ap of the basic fields and parameters F are found in the standard way (see e.g. 
[T6l [T7] ) and especially J27] for the analogous scalar problem) and have the forms 

A[v n ] = nA v = n(l -y/3), A[v']=d-A v , 
A[6 n ] = nA e = n(-l + y/6), A[9'} = d-A e , 

A U) = 2-y/3, A m = l. (2.12) 

These results are exact due to the exact expressions j u (g*, u*) = j K (g*, u*) = y/3 which 
follow from the relations (I2.8P at the fixed point. 

3. Inertial-range behaviour of the structure functions and the OPE 

Solution of the RG equations gives the asymptotic expressions for the various Green 
functions in the IR range, that is, for Ar ^> 1 and any fixed value of mr, where A is the 
UV momentum scale from ( 12. 31) and m is the IR scale from the correlators ( 11. 41) and 
( II. 7p . In particular, for the structure functions ( 11.11) one obtains: 

S 2n (r) = D^r- 2nA °Umr), (3.1) 
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cf. [27] for the scalar case. The inertial range corresponds to the additional condition 
mr <C 1. The asymptotic form of the scaling function ((mr) for small mr is determined 
by OPE and has the form 



where A^ are the critical dimensions of the relevant composite fields (composite 
operators) F and C n p(mr) are coefficients regular in (mr) 2 . 

Obviously, the leading terms of the asymptotic behaviour of the function ( 13. 21) for 
mr <C 1 are determined by the operators with smallest dimensions A^. However, some 
additional considerations should be taken into account. The operators whose dimensions 
appear in (13. 2 j) are those allowed by the symmetry of the model and by the symmetry of 
the quantity on left-hand side. In our model, these are scalar operators, invariant with 
respect to Galilean transformation and with respect to the shift 9 — > 9 + const. The 
number of the fields 9 in the operator F cannot exceed their number on the left-hand 
side: this is a consequence of the linearity of the original stochastic equation in 9. The 
operators which can be represented as total derivatives, dF, have vanishing mean values 
and do not contribute to (13. 2p . For more detailed discussion of all these points see e.g. 



In models of critical behaviour (like e.g. the 4 model) the leading contribution to 
the OPE is given by the simplest operator F = 1 with A^ = 0. The feature specific 
of the models of turbulence is the existence of the so-called "dangerous" operators with 
negative critical dimensions A^ < 0. They dominate the small- mr asymptotic behaviour 
of the scaling functions and lead to singular dependence of the IR scale, that is, to the 
anomalous scaling (51 [6] . 

Like in the scalar [27] and rapid-change vector [10] cases, the most dangerous in 
our model are the simple operators 9 whose dimensions are known exactly, see (I2.12p . 
But they are not invariant with respect to the shift and do not contribute to (13. 2 j) . Thus 
the leading terms for mr <C 1 in (I3.2p are given by the operators with minimal possible 
number of spatial derivatives, which guarantee the invariance with respect to that shift, 
that is, one derivative per each field. They have the forms 89 . . .89 with even number 
of factors 89 (otherwise it is impossible to get a scalar of the needed form; see below) 
and all possible contractions of the vector indices of the fields and derivatives. For a 
scalar field 9 there is only one variant of the contraction: (8i98i9) n , and A ra in f 1 1 . 2 1) is 
the critical dimension of this operator [3 |27J . 

For the vector field 9 the number of possible contraction variants in the structure 
(89) 2n rapidly increases with n. All the operators with a given n mix heavily in 
renormalization, so that the leading exponents A^ in H3 .21) are not determined by an 
individual operator, but rather are the eigenvalues of the matrix of critical dimensions. 
The minimal eigenvalue gives the leading term of the small-mr behaviour, and the others 
give the corrections which, for small y, can be very close to the leading term. 

For n = 1 the contraction variant is still unique: 8i9jdi9j. Like in the scalar case 
[5], its dimension A n = is found exactly from a certain Schwinger equation, so that 




(3.2) 



F 



[27\. 
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the function S2 oc r _2Afl reveals no anomalous scaling; cf. [TU] for the rapid-change case. 
(The second variant didjdjOi = didj(6j8i) leads to a total derivative and thus gives no 
contribution to (13. 21) ). However, for n = 2 there are already 7 variants [H HDJ, arid for 
n = 9 as many as 47 246 [T2] . 

Furthermore, there are some nontrivial linear relations between the operators with 
the same n, which reduce the number of independent monomials, especially for low 
values of d (up to 6 for n = 2 and general d and up to 154 for n = 9 and d — 3), but give 
rise to a difficult problem of finding and exclusion of all the redundant operators; see 
[T2] for a detailed analysis. In two dimensions, the transverse vector field is expressed in 
terms of a scalar field by means of the antisymmetric Levy-Civita tensor: 9{ = e^dkip- 
This allows one to diagonalize the matrix of critical dimensions of the operators (d6) 2n 
for arbitrary n and to derive the explicit results for the leading anomalous exponents to 
the order £ 2 [12]. 

For a general d, the families with different n should be studied separately. 
Surprisingly enough, the problem drastically simplifies for d — >■ 00, which allows one 
to achieve very high values of n and to obtain simple explicit results for leading and 
correction exponents [ID] 111]. So far, all these results were confined with the Kraichnan's 
velocity ensemble. 

4. Critical dimensions of the operators (d9) 2n for large d. 

The analysis of the renormalization of the composite operators (d8) 2n and the 
practical first-order (one-loop) calculation of the renormalization constants and critical 
dimensions in the model (12. II) . (12. 2 1) is very similar to the case of vector Kraichnan's 
model, which is discussed in great details in ref. [ID] ; see especially sec. B and C 
and appendix B. The only relevant one-loop Feynman diagrams in the two models 
differ only in the scalar factor stemming from the integrals over the frequency; all the 
tensor factors (projectors, vertices etc.) are exactly the same. As a consequence, the 
expressions for the renormalization constants here can be obtained by the substitution 
gSd/s — > gSd/yu{u + 1) in expressions like (6.19), (B15) and (B17) in [ID] , while the 
final fixed-point expressions for the critical dimensions are obtained by the replacement 
e — > y/3 in expressions like (6.24), (6.25) or (B6c) (e = £ in the notation of [10J). 

For this reason, below we only briefly discuss the renormalization of the operators 
(d6) 2n and the simplifications that occur for large d. Detailed justification given in 
ref. [10] for the rapid-change vector model literally applies to the present case. 

The critical dimension A of an arbitrary scalar composite operator of the form 
(d6) 2n in the first order of the expansion in y has the form A = A\(d)y + 0(y 2 ) with 
a certain coefficient Ai(d) that depends on d. The first terms of its expansion in 1/d 
have the forms 

A x (d) = 2k + A u /d + 0(l/d 2 ), 

where k is an integer number satisfying the inequalities < k < n and A n is a numerical 
coefficient independent of y and d. It turns out that in the first order in 1/d the subsets 
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with different k "decouples" from one another, so that their renormalization can be 
studied separately. 

It is clear that for large d, dangerous operators with An < can be present only 
in the subsets with k = 0. They are formed by the operators (d6) 2n of a very special 
type, namely, those in which all the fields are contracted only with the fields, and the 
derivatives are contracted only with derivatives. For a given n all such operators are 
represented as the products 

F=(<t> 1 r(for*...(<t> q ) n *, (4.i) 

where Y2k=i ^ n k — n and 0^ is a scalar operator that includes 2k factors 89 and cannot 
be reduced to a product of certain scalar factors. Such a basic factor can be written in 
the form 

fc = d h e Sk d h e 31 d h e sl d h e S2 d h e S2 d h e sz ■ ■ ■ ^"e 8h _^"e 8k . (4.2) 

Let us give a few examples. For n = 2 there are two operators of the type (14.11) : 

F = {<Pl <M, 

for n = 3 there are three operators: 

F = {0?, 0102, 3 }, 

for n = 4 there are five operators: 

F = {0i, 0102, 02, 0103, 04}, 

and for n between 5 and 11 the number of relevant operators equals to 7, 11, 15, 22, 30, 
42 and 56, respectively. These number are much smaller than the total numbers of the 
operators (d6) 2n with a given n: for example, 2 rather than 6 for n = 2 and 30 rather 
than 47 246 for n — 9. However, for n = 28 there are as many as 3718 operators of the 
type (14.11) . so that the problem remains highly nontrivial even for d — > oo. Let us give 
the whole set of the operators of the type (14.11) for n = 7: 

F = {(Pi <^02 , 0102 , 0102 , 0103 , 0103 , 0?04, 0?05, 

0106, 0?0203, 0203, 010204, 0205, 0304, fa} (4.3) 

and for n = 8: 

F = {0?, 0102, 0^, 0?02, 02, 0^03, 0?03, 
0203 , 0?0203, 010203, 04, 0104, 0?0204, 0204, 

010304, 0105 , 010205 , 0305 , 0?06, 0206, 0107, 0s}- (4-4) 

Renormalization of the families of operators of the type (14.11) with different n can 
be studied separately: due to the linearity of the original problem (11.31) . the operators 
(d9) 2n do not admix in renormalization to the operators (d8) 2k if n > k. The leading 
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term of the double expansion in y and 1/d for the matrix of critical dimensions of the 
family (14. ip with a certain given n has the form 

A— + (4.5) 

where the ellipsis stands for the corrections in y and 1/d and A is a matrix with non- 
negative integer elements, which are determined by certain simple rules [TO]. 

The diagonal element A aa , corresponding to a certain operator F a with a given n, 
is given by the expression 

q 

A aa = n - n x - ^ n k k(k - 1), (4.6) 

k=2 

where nk is the number of factors (pk that constitute the operator F a . 

The non-diagonal elements are determined by the "fusion and decay processes" of 
the simple factors 0fc. Let us choose in the operator F a a pair of the simple factors 4>k 
and 4> p with certain k and p (k = p is allowed) and replace them by a single factor <pk+p- 
Then the original operator F a turns to some other Fp with the same n. This "fusion 
process" 0fc0 p ~^ <t>k+p gives to the matrix element A a p a contribution of the form Akp 
with the summation over all possible pairs of factors 4>k4>p entering into F a . For example, 
starting with the operator F± = <p\ one can get the operator F 2 = 0102 by means of the 
fusion 0i0i — > 02 with k = p = 1; this gives the matrix element A 12 = 4kpC^ = 12, 
where the factor C\ = 3 arises as the number of possibilities to choose the pair 0i0i 
from the three factors of the type 0i in Fx. 

Furthermore, let us choose in the operator F a some simple factor 0^ and replace it 
with the factor (pk- P 4>p with a certain 1 < p < (k — 1). Then the new operator Fp with 
the same n results. This "decay process" 0^ — > 0fc_ p p gives to the matrix element A a p 
a contribution of the form 2k per each factor 0^ entering into F a , if the factors 0fc_ p and 
P are different, that is, p ^ k — p, or k per each factor 0^ if they are identical, that is, 
p = k — p. 

Example: two possible decays in the operator F\ = 04 {k = 4) give rise to the 
operators F 2 = $\ [k = p — k = 2) and F 3 = 0103 (p = 1, p — k = 3); the corresponding 
matrix elements equal to A12 = k = 4 and A13 = 2k = 8. Another example: the decay 
02 — > 0i0i in the operator Fi = 0| gives rise to F 2 = 0i0| with the matrix element 
A12 = 3k = 6; the factor 3 accounts for the presence of three monomials 02 in the initial 
operator 

If the operator F a gives rise to another operator F@ as a result of certain fusion 
0fc0?> — > 4>k+ P , then it is clear that Fp gives rise to F a as a result of the "inverse decay" 
4>k+ P — >• 0fc0 P - Thus the matrix elements A a p and Ap a can vanish only simultaneously 
(which happens very often); otherwise they both are not equal to zero (but, in general, 
are not equal to each other). 
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2 



For n = 2 and 3 the matrices A have the forms: 

12 0\ 

2 8 1. (4.7) 
.0 6 3/ 

For the families of operators that mix in renormalization the exponents A^ in (13. 2 j) 
are determined by the eigenvalues of the matrices A with the subsequent substitution 
into (14. 5ft . The eigenvalues of the matrices (14. 7ft for n = 2 are ±2^2, while for 
n = 3 they are equal to 1 + 10 cosip = 9,673557, 1 — 5 cos-?/' — 5\/3 sinip = —7,64689, 
1 — 5 cos ip + Sv^sin-^ = 0,973333, where we have denoted if) = (l/3)arctg (6\/ 434) . 
For higher n the eigenvalues were found numerically. 



The matrices A for n as high as 6 can be found in [10]. As two more examples, 



;ive 


the matrices for the 


case n 


= 7: 
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14 


35 J 



and for n = 8: 
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In the both cases the operators are numbered according to their order in ( 14. 3 p and 

gap. 

The algorithm described above for constructing the matrices A was realized as a 
computer program; using it we have found all the matrices and their eigenvalues up 
to the family with n = 28 (it involves 3718 relevant operators). The computations for 
larger n become too time compiling (mostly for finding the eigenvalues) but they appear 
unnecessary: the eigenvalues demonstrate interesting regularities, which allows one to 
suggest for them some simple "empiric" formulae that become nearly exact for n large 
enough. 

5. Critical dimensions of the operators (d9) 2n for large n. 

The analysis of the matrices A with n < 28 reveals the following properties. 

All the matrices can be brought to diagonal form; all their eigenvalues are real and 
differ from zero. (One can check that the low-order matrices A can be made symmetric 
by a proper normalization of the basis operators (14.11) ; this is likely true for all orders. 
However, the elements of the matrices become non-integer after that procedure.) Among 
the eigenvalues for a given n there are always positive and negative. Starting from n = 3, 
the number of positive eigenvalues exceeds the number of negative eigenvalues (roughly 
speaking, twice). The maximum (by modulus) positive and negative eigenvalues grow 
monotonously with n, and the maximal positive eigenvalue is always larger (by modulus) 
than the negative one (roughly speaking, twice for large n). Let us give a few examples. 
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1575; 


1080.5; 


-574.12 


for 


n 


= 24, 


1958; 


1175.5; 


-623.11 


for 


n 


= 25, 


2436; 


1274.5; 


-674.11 


for 


n 


= 26, 


3010; 


1377.5; 


-727.11 


for 


n 


= 27. 



For n = 28 there are 3718 operators of the form fl4.ll) . 2569 positive eigenvalues and 
1149 negative ones. The maximal positive eigenvalue of the matrix A equals to 1484.5, 
the minimal one is —782.1. For smaller values of n the same five numbers are: 

1072; 503; 

1337; 621; 

1674; 762; 

2070; 940; 3010; 1377.5; -727.11 forn = 27. (5.1) 

The relation (14. 5 p shows that dangerous operators with negative critical dimensions 
correspond to positive eigenvalues of the matrices A, which therefore are the most 
interesting ones. Furthermore, the maximal (for a given n) eigenvalue determines the 
leading term in the IR asymptotic behaviour of the structure function S 2n in ( 13. ip , ( 13. 2 p , 
that is, the principal anomalous exponent. The other eigenvalues determine corrections 
to the leading term; they diverge for mr — > if the corresponding eigenvalue is positive 
and decrease if it is negative. 

The operators that possess definite critical dimensions ("scaling operators") are 
certain linear combinations of the basis monomials ( 14. ip . It turns out that the scaling 
operators corresponding to the maximal eigenvalues involve all the monomials from 
the family with the given n; all the coefficients are positive and become closer to each 
other as n grows. (All the other scaling operators necessarily involve coefficients with 
different signs, because the eigenvectors of the matrix A in its symmetric form must be 
orthogonal.) This "democracy of monomials" can be opposed to the two-dimensional 
case, where the principal eigenvalues correspond to the operators of a very special form: 
powers of the local dissipation rate (c^c^ — c^c^j); see [12]. 

Let us give all maximal positive eigenvalues Xo{n) of the matrices A for all n from 
2 to 28: 

2.828; 9.67356; 20.617; 35.5888; 54.5717; 77.5602; 

104.5518; 135.55; 170.54059; 209.5366; 252.5334; 

299.53063; 350.52832; 405.53; 464.5246; 527.52308; 

594.52175; 665.52055; 740.51949; 819.51852; 902.51765; 

989.51686; 1080.5; 1175.5; 1274.5; 1377.5; 1484.5. (5.2) 

The figure 1 shows that the eigenvalues Ao(n), plotted against the number n, are 
approximated nicely by a smooth curve (the upper solid line). Surprisingly enough, that 
curve is described very well by a simple analytic expression: 

A (n) = 2n 2 -3n + 1/2 + 0(l/n), (5.3) 

where the 0(l/n) correction is rather small already for modest values of n. Indeed, the 
inspection of the maximal eigenvalues Xo(n) given in (15. 2p shows that the expression 
2n 2 - 3n gives exactly (and with no exceptions) their integer parts, and the refined 




Figure 1. The eigenvalues of the matrices A. The circles correspond to the leading 
branch Ao from (|5.3[) . the squares, triangles and asterisks correspond to the branches 
Ai, A2* and A3, from (|5.4I) . 



expression 2n 2 — 3n + 1/2 gives, starting from n = 5, the first number after the decimal 
point (it is equal to 5 for all n > 5). 

The next-to-maximal eigenvalues also form a smooth curve, the next ones form 
their own branch, etc. All the branches are described by simple explicit formulae, which 
rapidly become nearly exact as n grows. Let us give them for a few branches closest to 



the leading one (j5.3j) 






\ (n) 


= 2n 2 


- 3n + 1/2, 


Ai(n) 


= 2n 2 


-7n + 7/2, 


A2,l( n ) 


= 2n 2 


- 11/2 + 40/3, 




= 2n 2 


- lln + 23/3, 


A3,l( n ) 


= 2n 2 


- 15n + 187/6, 


A 3 ,2(» 


= 2n 2 


- 15n + 135/6, 


A 3 ,3(n) 


= 2n 2 


- 15n + 83/6, 


A 4 ,i(n) 


= 2n 2 


- 19n + 57.1, 


A 4) 2(n) 


= 2n 2 


- 19n + 44.3, 



with corrections of order 0(l/n). 
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It is interesting to note that, for n = 1, the only eigenvalue A(l) = (which is known 
exactly) belongs to neither curve. Probably this fact is related to the observation that 
the agreement between the formulae (15. 3ft . (15 Ah and the exact numerical values for the 
eigenvalues improves if the correction term is written in the form 0(l/(n — 1)). 

In their turn, the expressions ( 15 .31) . H5.4I) demonstrate interesting regularities: the 
leading (quadratic in n) term is the same for all the branches; the next-to leading term 
(linear in n) is negative with the odd coefficient growing with the step of 4. The branches 
can be grouped according the value of that coefficient, and the number of branches with 
a given value grows: there is one branch in the first two groups, two branches in the 
third group, three branches in the fourth group etc. It is also worth noting that if the 
constant (independent of n) terms in (I5.3P are replaced with the closest integer numbers 
and the 0{l/n) corrections are neglected, the resulting expressions give exactly the 
integer parts for all the eigenvalues. 

These regularities are also illustrated by figure 1, where all the positive eigenvalues 
of the matrices A are shown with n from 3 to 15; according to (14. 5ft . they correspond 
to "dangerous" composite operators with negative critical dimensions. The solid lines 
correspond to representatives of the principal branches: Ao from (15. 3j) and Ai, A24, 
A 3>3 from (15.41) . They are plotted according the formulae (15.31) and (15.41) neglecting the 
0(l/n) corrections. The other branches from the groups A2,* and A3.* are not shown by 
solid lines in order to make the picture more graspable. The circles denote the maximal 
eigenvalues of the matrices A, and the squares denote the next-to-maximal ones; they 
lie exactly on the principal branches Ao and Ai from (15. 3ft . (15.41) . The eigenvalues that 
correspond to the two branches of the group A 2 ,* are denoted by triangles and the 
eigenvalues of the three branches of the next group A3,* are denoted by asterisks. 

The general formulae (\5.3h , (15.41) for negative eigenvalues do not work so well in 
comparison with positive ones. This fact is illustrated by the same figure 1, where some 
negative eigenvalues (with n from 4 to 6) are also shown. It turns out, however, that the 
negative eigenvalues form their own pronounced branches; the principal one is described 
by the empiric formula 



Comparison of expressions (15.31) and (I5.5P shows that the ratio of the maximal positive 
and maximal (by modulus) negative eigenvalues of the matrix A tends to 2 as n grows, 
in agreement with the numerical values for 24 < n < 28 given in (15.21) . 

6. Conclusion 

By means of the field theoretic renormalization group and operator product expansion 
we studied the problem of anomalous scaling of the transverse vector field passively 
advected by a turbulent velocity field. The dynamics of the vector field is governed by 
the stochastic equation (11.31) . (11. 4ft . while the velocity was described by the stirred 
NS equation (11.61) . (11.71) . (11.81) . The anomalous scaling arises as a consequence of 




n 2 + 2 + 0(l/ra). 



(5.5) 
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the existence in the OPE of the so-called "dangerous" composite fields (operators) 
with negative critical dimensions. The leading terms of the inertial-range asymptotic 
behaviour of the structure functions ( 13. lj) . (13.21) are determined by the matrices of 
critical dimensions of the families of composite fields of the form (d8d8) 2n . For d — > oo, 
dangerous operators can be present only in the subsets of operators of the special form, 
(I4.ip . (14. 2p . and the corresponding matrices of critical dimensions can be constructed 
by means of a simple algorithm. This allowed us to calculate them in the leading order 
of the double expansion in y and 1/d up to the order n = 28. The eigenvalues of those 
matrices (that is, the critical dimensions of the corresponding families of operators) 
demonstrate intriguing pronounced regularities. This fact allows one to describe them 
by simple empiric formulae which become practically exact as n grows. In particular, the 
leading term of the asymptotic behaviour of the structure function ( 13. If) in the inertial 
range has the form 



there are also explicit expressions for the correction exponents. Thus, the complete 
description of the anomalous scaling for our vector model is given for all n. 

The regularities demonstrated by the critical dimensions of the sets of composite 
fields and their branches, discussed in the present paper, are so intriguing that we 
cannot but think that some unknown symmetry lies behind them. In this connection 
it is interesting to note that the critical dimensions of certain composite operators in 
quantum chromodynamics (QCD) and in the N = 4 supersymmetric gauge theories also 
show interesting behaviour (also in the large- n limit and in the one-loop approximation): 
in particular, for the QCD case the corresponding evolution equations appear equivalent 
to the integrable Heisenberg model [28] . 

The Kraichnan's rapid-change model of passive scalar advection is sometimes 
referred to as "the Ising model of fluid turbulence." In this connection, it is worth 
recalling that the original Ising (or, better to say, Lentz-Ising) model of magnetism, first 
introduced in 1920 [29J, has still remained a source of inspiration for new physical and 
mathematical ideas and techniques like integrability, fermion-boson transformations, 
conformal invariance and discrete holomorphisity: for a recent discussion, see [30j and 
references therein. 

One may think that, in spite of a great deal of work devoted to Kraichnan's model 
and its descendants, the deep physical and mathematical contents that lie behind them 
are not completely disclosed. We believe that identifying the hypothetical symmetry of 
the passive vector problem that gives rise to the intriguing regularities discussed in the 
present paper will help to reach a deeper understanding of the anomalous scaling in the 
real fluid turbulence. 
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